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I. INTRODUCTION 



A long standing question for the embedded ensembles is about their analytical tractability. 
Amenability to mathematical treatment is one of the four conditions laid down by Dyson [1] 
for the validity of a random matrix ensemble. Simplest of the two-body unitary ensemble is 
the embedded Gaussian unitary ensemble of two-body interactions [EGUE(2)] for spinless 
fermion systems. For m fermions in N sp states, the embedding is generated by the SU (N) 
algebra and although this ensemble is known for many years, only recently [2], after the 
first indications implicit in [3, 4], it is established that the SU(N) Wigner-Racah algebra 
solves EGUE(2) and also the more general EGUE(fc) [as well as EGOE(fc)]. These results, 
with U(N) algebra, extended to BEGUE(A;) for spinless bosons in N sp states (see [2, 5]). 
For EGUE(2)-s for fermions with spin and EGUE(2)-S77(4) for fermions with Wigner's 
spin-isospin SU(4) symmetry, the embedding algebras, with Q number of spatial degrees 
of freedom for a single fermion, are U(Q) <8> SU(2) and U(Q) <S> SU(4) respectively. It 
was shown in [6, 7] that the Wigner-Racah algebra of these embedding algebras will allow 
one to obtain analytical results for the lower order moments of the one- and two-point 
correlation functions in eigenvalues. Similarly, following the recent work [8, 9] on BEGOEs, 
it is easy to recognize that the embedding algebras for BEGUE(2)-F for two-species boson 
systems with F-spin and BEGUE(2)-S77 (3) for spin one boson systems are U (Q)<S>SU (2) and 
U (Q)®SU (3) respectively. The purpose of the present paper is to establish on one hand that 
the Wigner-Racah algebra of these embedding algebras solve the corresponding embedded 
unitary ensembles and on the other to generalize the formalism to any EGUE(2) with U (Q)<g> 
SU(r) embedding and generated by random two-body interaction with SU(r) symmetry. 
Hereafter we call these ensembles EGUE(2)-S77(r) and they apply to both fermion and 
boson systems. 

In Section 2, given is the general formulation based on Wigner-Racah algebra for lower 
order moments of the one- and two-point functions in eigenvalues generated by EGUE(2)- 
SU(r) (r is any positive integer, r > 1). Sections 3, 4 and 5 give analytical results for boson 
systems with r = 1, r = 2 and r = 3 respectively. In addition, some numerical results for 
lower order correlations generated by these ensembles are also given in Section 5. Finally, 
Section 6 gives concluding remarks. 
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II. EGUE(2)-SE/(r) ENSEMBLES: GENERAL FORMULATION 



Consider a system of m fermions or bosons in f2 number of sp levels each r-fold degenerate. 
Then the SGA is U(rQ) and it is possible to consider U(rQ) D U(Q) <g> SU(r) algebra. 
Now, for random two-body Hamiltonians preserving SU(r) symmetry, one can introduce 
embedded GUE with U(Q)<g>SU(r) embedding and this ensemble is called EGUE(2)-S77(r). 
Ensembles with r = 2, 4 for fermions correspond to fermions with spin and spin-isospin 577(4) 
symmetry. Similarly, for bosons r = 2, 3 are of interest. Also r = 1 gives back EGUE(2) and 
BEGUE(2) both. It is important to note that the distinction between fermions and bosons 
is in the U (Q) irreps that need to be considered. Now we will give a formulation in terms of 
SU(£l) Wigner-Racah algebra (the SU(r) algebra involved will be simple as H has SU(r) 
symmetry) that is valid for any r. The discussion in the remaining part of this Section is 
essentially from [7] but it is repeated briefly not only for completeness but also to generalize 
it to any r and also to bosons systems (in [7], fermions with r = 4 is used). 

Let us begin with normalized two-particle states \f2F2', V2P2) where the U(r) irreps 
F2 = {l 2 } and {2} and the corresponding U(Q) irreps f 2 are {2} (symmetric) and {l 2 } 
(antisymmetric) respectively for fermions and {l 2 } (antisymmetric) and {2} (symmetric) 
respectively for bosons. Similarly V2 are additional quantum numbers that belong to f 2 
and 02 belong to F 2 . As f 2 uniquely defines F 2 , from now on we will drop F 2 unless it 
is explicitly needed and also we will use the f 2 -H- F 2 equivalence whenever needed. With 
A^(f 2 v 2 (3 2 ) and A(f 2 v 2 (3 2 ) denoting creation and annihilation operators for the normalized 
two particle states, a general two-body Hamiltonian operator H preserving SU (r) symmetry 
can be written as 

H = H {2} + 77 {12} = H hvH {2)A\f 2 v{h)AU2V^ 2 ). (1) 

/ 2 ,*4,4,/3 2 ;/ 2 ={2},{l 2 } 

In Eq. (1), H^ v i v f(2) = (^f 2 v(f3 2 \ H \ f2V t 2 (3 2 j independent of the (3 2 s. The uniform sum- 
mation over 02 in Eq. (1) ensures that H is SU(r) scalar and therefore it will not connect 
states with different /2's. However, H is not a SU(r) invariant operator. Just as the two 
particle states, we can denote the m particle states by |/m^m^m)i F m = f m for fermions and 
F m = fm for bosons. Action of H on these states generates states that are degenerate with 
respect to f3^ but not v f m . Therefore for a given / m , there will be dn(f m ) number of levels 
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each with d r (f m ) number of degenerate states. Formula for the dimension dn(f m ) is [10], 

n 



Mfm) = n 

i<j=l 



fi-fj + j 

j-i 



(2) 



where, f m = {/i,/2, •••}■ Equation (2) also gives d r (F m ) with the product ranging from 
i — 1 to r and replacing /, by Fj. As if is a SU (r) scalar, the m particle H matrix will be a 
direct sum of matrices with each of them labeled by the / m 's with dimension dn(f m ). Thus 



fm 



(3) 



It should be noted that the matrix elements of Hf m (m) matrices receive contributions from 
both H {2} (2) and H {l2y (2). 

Embedded random matrix ensemble EGUE(2)-S77(r) for a m fermion or boson system 
with a fixed f m , i.e. {Hf m (m)}, is generated by the ensemble of H operators given in Eq. 
(1) with iJ{ 2 }(2) and ii{i2}(2) matrices replaced by independent GUE ensembles of random 
matrices, 

{if (2)} = {H {2} (2)} GUE © {H {1 2 } (2)} GUE . (4) 

In Eq. (4), { } denotes ensemble. Random variables defining the real and imaginary 

parts of the matrix elements of if/ 2 (2) are independent Gaussian variables with zero center 
and variance given by (with bar representing ensemble average), 



H f2v\vl^) H f ! 2V 3 v i(2) = 5f 2 f^5 v l v 45 v 2 v 3 (A /2 ) . 

Also, the independence of the {if{ 2 }(2)} and {if{i2}(2)} GUE ensembles imply, 



(5) 



#{2KS 2 ( 2 ) 



H {1 2 }V 3 V 4(2) 



Q 



#{2K% 2 (2) 



H {1 2 }V 3 V ^(2) 



Q 



(6) 



for P and Q even and zero otherwise. Action of H defined by Eq. (1) on m particle 
basis states with a fixed f m , along with Eqs. (5)- (6) generates EGUE(2)-5'?7(r) ensemble 
{Hf m (m)}; it is labeled by the U(Q) irrep f m with matrix dimension dn(f m ). 

As shown in [2, 6, 7], tensorial decomposition of H with respect to the embedding algebra 
U(Cl)®SU(r) plays a crucial role in generating analytical results; as before U (Q) and SU(Q) 
are used interchangeably. As H preserves SU(r), it transforms as the irrep {0} with respect 
to the SU(r) algebra. However with respect to SU(Q), the tensorial characters, in Young 



su(n) 



[2] 



{2"-'} 



{2"} {32«- 2 l} {42 n - 2 } 

= {0} = {21<«} 

v=0 v =l v=2 



(a) 



SU(Q) 



1 ' 2 ! 



{!<«} 



£2-4 



{21<«} 

= {0} 

v=0 v=l 



{2 2 1«~ 4 } 
v=2 



(b) 



FIG. 1. Young tableaux for the tensorial parts of a two-body Hamiltonian with respect to SU(£l) 
algebra. Young tableaux for various (a) tensorial parts with respect to SU(£l) for the f 2 = {2} 
part of H; (b) tensorial parts with respect to SU(Q.) for the f 2 = {l 2 } part of H. 

tableaux notation, for f 2 = {2} are F v = {0}, {21 n ~ 2 } and {42 n ~ 2 } with v = 0, 1 and 
2 respectively. Similarly for f 2 = {l 2 } they are F u = {0}, {21 n " 2 } and {2 2 1^- 4 } with 
v — 0,1,2 respectively. Note that F v = f 2 x / 2 where / 2 is the irrep conjugate to / 2 and 
the x denotes Kronecker product. Given a U(Q) irrep {/} = {fi, f 2 , ■ ■ ■ , fn}, we have 
/ = {/i — /n, /i — /n-i, • • • , fi — f 2 , 0}. Young tableaux for the F„'s are shown in Fig. 1. 
Now, we can define unitary tensors 5's that are scalars in SU(r) space, 

B(f 2 F u u u ) = A\f 2 v f 2 (3 2 )A(f 2 v l 2 (3 2 ) (f 2 v f 2 J 2 4\ F u u u ) 

vi4,p2 (7) 

x (F 2 (3 2 T 2 J 2 | 00) . 

In Eq. (7), (f 2 ) are SU(Q) Wigner coefficients and {F 2 ) are SU(r) Wigner 

coefficients. The expansion of H in terms of 5's is, 

H= Yl W{f 2 F v u v )B{f 2 F v u v ) . (8) 
/2,-FV, u v 

The expansion coefficients W's follow from the orthogonality of the tensors -B's with respect 
to the traces over fixed f 2 spaces. Then we have the most important relation needed for all 
the results given ahead, 

WU^WU^K) = S hfa 6 F F <W (X h ) 2 d r (F 2 ) . (9) 
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This is derived starting with Eq. (8) and using Eqs. (4)- (7) along with the sum rules for 
Wigner coefficients appearing in Eq. (7). 

Turning to m particle H matrix elements, first we denote the U(Q) and U(r) irreps by 
f m and F m respectively. Correlations generated by EGUE(2)-5'?7(r) between states with 
(to, f m ) and (m',f m >) follow from the covariance between the m-particle matrix elements 
of H. Now using Eqs. (8) and (9) along with the Wigner-Eckart theorem applied using 
SU(Q) <S> SU(r) Wigner- Racah algebra (see for example [11]) will give 



f m F m vLP | H I fmFmVl/3) UmiFm'V S m ,V \ H \ fm'F m ^,(5> 



^2 i ff) ^ ^ m "' B ^ 2F ^ "' "' B (h Fv ) III f m '} P ' 

/ 2 ,.F v ,u v P'P 



\2 



X (fmVin F V U V | f m vQ o (fm'V^ F v bJ v \ f m ,V f m , / ^ , 



(fm HI B{f 2 F v ) HI f m ) = y F(m) U (fmhfmf2- t fm-2F v ) p 

Here the summation in the last equality is over the multiplicity index p and this arises as 
f m x F v gives in general more than once the irrep f m . In Eq. (10), 

F(m) = -m{m- l)/2 , (11) 

dn(fm) is given by Eq. (2) and (. . .) and £/(•••) are SU(£l) Wigner and Racah coefficients 
respectively. Similarly, A// m is dimension with respect to the S m group [10], 



p 

ml 



n & - « 



^>-= ^i-v 1 ^a*'-'- < 12) 

Note that p denotes total number of rows in the Young tableaux for f m . 

Lower order cross correlations between states with different (to, f m ) are given by the 
normalized bivariate moments T> P q (to, f m : to', f m '), P — Q — 1, 2 of the two-point function 
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S p where, with p m ^ m (E) defining fixed-(m, f m ) density of states, 



S mfm:m ' f ™'(E,E') = p m 'fm(E)p m ''fm'(E') - p m >f™(E) p m '>f m >(E') ; 

E u (m, f m : m', / m ,) = (H) mJ ™ (H) m '^'/^ (H*)^ (tff^ , (13) 



\m,f m /TT2\ m '<f n 



- 1 . 



In Eq. (13), (H 2 ) mJm is the second moment (or variance) of the eigen value density p m >f™(E) 
and its centroid (H) m ' fm = by definition. As (H) mJm is the trace of H (divided by 
dimensionality) in (m, f m ) space, only F v = {0} will generate (H) m ' fm (H) m '' s ™'. Then 
trivially, 



{H) m,f m {H) m>J m , = ^^p/ 2 ( m j m )p/ 2(ffl 'j m ,) ; 

72 



(14) 



fm-2 



Mr 

Jm 



Writing ^Jj 2 ) m ^ m explicitly in terms of m particle H matrix elements, 



and applying Eq. (10) and the orthonormal properties of the SU(Vt) Wigner coefficients 
lead to 

^ W?) ^ 2 

where 



Q u (f2--m,f m ) = [F(m)] 2 E !r 2 x uu(f2\f m -2,fL-2,*'v)- (16) 

Jm-2;J m _2 



The Xj/tr function involves S'C/ (f2) Racah coefficients, 

Xuu(f2', fm-2, f' m -2, & v) = 

y-v (f m , f 2 , fm, fm-2, F ' u ) p U (f m , f 2 , f m , /gj /m-2; -^V)p (1?) 
p U (f m , f 2 , fm, f2, fm-2, {®})U(f m , f 2 , fm, f 2 \ f' m - 2 , {0}) 

Summation over the multiplicity index p in Eq. (17) arises naturally in applications to 
physical problems as all the physically relevant results should be independent of p which is 
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a label for equivalent SU(Q) irreps. It is easy to see that, 

Q u=0 (f2:m,f m )= [Pf\mj m )] 2 . (18) 

Eqs. (14)-(16) and Table 4 of [7] will allow us to calculate covariances En in energy centroids; 
Table 4 of [7] is a simplified version of the tables in [12]. For the covariances £ 22 in spectral 
variances, the formula is [7] 

V (m f .„/ f \ _ X W + X (i 2 } + 4X {i 2 l{2} . 
2->22{m, f m , m , j m ,) - — — , 

X ^ = TJTT^ 2 E [^(^r 1 fi"(/2 : ™>f m )&(f2 "'Jm') , (19) 
lrfn(/2jj ,,=0,1,2 

X{12}{2} = d(^MOT) E [W] _1 »>,/™) k v WM ■ 

Here d^(F u ) is the dimension of the irrep .F„, and we have dn({0}) = 1, cfo({2, l n ~ 2 }) = 
Q 2 -l, d n ({4,2 n - 2 }) = fi 2 (fi + 3)(fi-l)/4, and^({2 2 ,l Q - 4 }) = Q 2 (Q - 3)(Q + l)/4. Note 
that Q I/ (f2 '■ m, fm) are defined in Eq. (16). The functions R u (m,f m ) also involve SU(Q) 
[/-coefficients, 

if(m,/ m ) = [F(m)] 2 £ ^- 2 ^- 2 y^(/^- 2 ,/;- 2 ;^); 

Jm-2,J m _2 

E g (/m, {1 ~ 2 }, /m, {I 2 }; Jm-2, -FY)p ^ (/m, {2"' 1 }, fm, {2}; /^_ 2 , 
W m Jl n - 2 }J.M{l 2 };/ m -2,{0})W m ,{2^},/ m ,{2};^_ 2 ,{0}) • 

In Y uu (f m _ 2 , f' m _ 2 ] F v ), f m -2 comes from f m <g> {l n ~ 2 } and /^_ 2 comes from f m <g> {2 n_1 }. 
Similarly, the summation is over z/ = and 1 only as v = 2 parts for / 2 = {2} and {l 2 } are 
different. Formulas for Yuy are given in Table 7 of [7] and they are simplified version of the 
results in [12]. It is useful to note that, 

7C=\m, f m ) = pW(m, fm)P {l2 \m, f m ) . (21) 

Compact analytical results collected in Tables 4 and 7 of [7] for Xjju an d Yi/f/ and Eqs. (2), 
(12) - ( 21) will allow one to derive analytical/numerical results for spectral variances and 
covariances in energy centroids and variances for any EGUE(2)-S , /7 (r) for fermion or boson 
systems. 
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Spinless bosons 



U(Q) irreps 



(i) 



fm 


={m} fm— 

m — 


2 ={m-2} 
















• 


• 



a a 



in 



(ii) 



Spinless fermions 

fm=a m } 

f m - 2 ={l m - 2 } 



FIG. 2. Young tableaux denoting the SU{VL) irreps f m = {m} and {l m } as appropriate for (i) 
spinless boson and (ii) spinless fermion systems. Removal of two boxes generating m — 2 particle 
irreps f m -2 for these systems are also shown in the figure. For (i) only the irrep f2 = {2} will 
apply and similarly for (ii) only {l 2 } will apply. 



III. RESULTS FOR BEGUE(2): r = 1 



Simplest of the EGUE(2)-5'f/(r) are the EGUEs with r = 1 and they corresponds to 
EGUE(2) and BEGUE(2) depending on totally antisymmetric or symmetric f m one consid- 
ers. Also they correspond to k = 2 in [2] and [5] for fermion and boson systems respectively. 
As detailed results for fermion systems are available in [2, 6, 8], in the present Section and 
in the next two Sections we consider only boson systems. Let us begin with BEGUE(2). 
For this ensemble, in order to apply the formulas given Section 2 for (H 2 ), £ n and S 22 , first 
we need formulas for Xuu and Yu\j. Some of these, taken from Tables 4 and 7 of [7], are 
given in Table I by reducing them to much small number of formulas. For applying these 
formulas, we need the 'axial distances' r^- for the boxes i and j in a given Young tableaux. 
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Given a f m = {/i, / 2 , . . . , /n} we have, 



r ij = fi-f j +j-i. (22) 



In terms of the functions Il£ 6 \ 11^ , n£ 6c \ 11^ and n" are defined as, 



ni 6) = n (! - V^) , 

i=l,2,...,Q;i^a,i^b 

ni a) = n ( x - ^ , 

i=l,2,...,f2;i^a,i^6 

ni 6c) = n C 1 - ; a^b^c, (23) 

i=l,2,...,Q;i^a,i^b,i^c 

k= n ( x - v^«) , 

i=l,2,...,f2;i^a 

= n c 1 - • 

i=l,2,...,f2;i^a 



With these we can calculate and 1W; see [7] for full discussion. For BEGUE(2), the 
algebra U(£l)®SU(r) with r = 1 reduces to just U(Q) or SU(Q). Similarly, / m is the totally 
symmetric irrep {m} and f m -2 = {m — 2}. Therefore to generate f m -2 only the action of 
removal of {2} from f m is allowed. Denoting the last two boxes of f m by a and a (note that 
we can remove only boxes from the right end to get a proper Young Tableaux and also boxes 
in a given row must have the same symbol to apply the results in Table I) as shown in Fig. 
2, we have 



IT' = 



m + i — 1 , 
m 



m + Q - 1 



(m + n- l)(m + n-2) 
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Similarly Aff m = 1 and A/ m _ 2 = 1 as both are symmetric irreps. Now the formulas in Table 
I will give Xjjjj and then using Eq. (16) we have, 

Q^({2};m,{m}) = m ^ m ~ lf , 

Q (I 2 )! ™> W) = , ^ol^ ' (25) 



Q- 2 ({2};m,{m}) = 



m(Q + 2) 

m 2 (m - l) 2 fi 2 (fi - 1)(Q + m)(Q + m + 1) 



4 2(Q + 2)m(m-l) 

These and Eq. (15) will give, 



<tf 2 > {m} = A ?2} (;)( fi + ;- 1 )^A 2 2} A-(^, m ,2); 



(26) 



A l/ (fi,m, A;) 



m — zA /T2 + m + z/ — 1 

fc J V fc 



This agrees with the result stated in [2, 5]. As P^(m,{m}) = —m(m — l)/2, we have 
easily, 

En({m},{m'}) = 2> /m(m - l)(m')(m' - 1) 

Q(Q + 1) v 7 ^ + ™ - + ™ - 2) (ft + m' - + m' - 2) 

Again this agrees, for m — m! with the result stated in [2, 5]. Further, S 2 2 is determined 
only by X{ 2 } defined in Eq. (19) and then, using Eq. (25), we have 

2 

£ 22 ({m}, {m 1 }) = — p — ■ 



36 l 3 iSl + 3 \ 2 /V 2 



X 



^9/^ x /fi + m+ 1\ /fi + m' + 1\ n9/ ^ Jm\fm'\ ( 28 ) 

4ff(«-l)( 2 j( 2 j + 4 (!i + 2«n + 3)( 2 j(^ 



+4(fi 2 - + 3)(m - + m)(m' - + m')] • 
For m = m', it can be verified that Eq. (28) reduces to 

and this agrees with the result given in [5]. Note that F u is {0}, {21 n_2 } and {42 n ~ 2 } for 
v — 0, 1 and 2 respectively. It is useful to mention that Eqs. (27) and (28) follow from the 
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results for fermion systems given in [13] with Q — > —Q symmetry. Finally, it is useful to 
mention that in the m — > oo and Q finite limit we have, 

2 

£n({m}, {m 1 }) 



(30) 

f ,r, r m n, o g!g -!) + (» + 2) 2 (^ + 3) + 4(^ 2 - l)(n + 3) 
£ 22 ({m}, {m }) = 8 fi2 (n + i)2 (n + 2) 2 (n + 3) • 

Non-vanishing of En and S22 for finite Q in the m — > 00 limit is interpreted in [5, 14] as 
non-ergodicity of BEGUE ensembles. See the discussion in [15] for the resolution of this 
problem. 



IV. EMBEDDED GAUSSIAN UNITARY ENSEMBLE FOR BOSONS WITH F- 
SPIN: BEGUE(2)-SU(2) WITH r = 2 

For two species boson systems we have BEGUE(2)-S'?7(2) and then the formulation in 
Section 2 with r = 2 will be applicable. Here the two species are assumed to be the two 
components of a fictitious F-spin as discussed recently in [8]. For such a m boson system, 
the SU (Vt) irreps will be two rowed denoted by f m = {m — r, r} with F = y — r. With 
this, there are three allowed f m ~2 irreps as shown in Fig. 3. The irreps in (i) and (iii) in the 
figure can be obtained by removing f 2 = {2} from f m . However for (ii) in the figure both 
{2} and {l 2 } will apply. For f m ~2 = { m — r — 2,r} irrep [this corresponds to (i) in Fig. 3] 
we have 

T a 2 = m - 2r + 1 , 

r ai = m — r + i — 1; % — 3, 4, . . . , , 

n , _ (m - 2r)(m - r + 1) (31) 



n'' 



(m — 2r + l)(m — r + f2 — 1) 

(m - 2r - \){m - r){m - r + 1) 



(m - 2r + l)(m - r + Q - l)(m - r + Q - 2) ' 
Similarly for j m -2 = {m — r, r — 2} irrep [this corresponds to (iii) in Fig. 3] we have 

7&1 = 2r — m — 1 , 

T~bi = r + % — 2 , % — 3, 4, . . . , Q 

(r)(2r-m-2) (32) 

6 ~ (2r-m- l)(r + fi-2) ' 
= (2r-m-3)(r)(r-l) 

6 (2r - m - l)(r + - 2)(r + - 3) ' 
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TABLE I. Formulas for X uu (f 2 ; f m -2, f' m -2\ F ») and Y Uu(fm-2, f m - 2 ; F v) witn v = 1,2. 



{fm-2} {f'm-2} 



X f7[/ ({l 2 };/ m _ 2 ,/^_ 2 ;{2-,l^}) 



{f(ab)}{f(ab)} 



g h 1 (»-!)(» -2) - 1) 



1 + 



1 



1 + d- 1 



Tab/ II 



r(«) 



1 4 A 



T ab J U ( a b) ^ 



{f(ab)}{f(ac)} 



0(0-1) 2 . 4 r , n , 1 



{fm-2} {fm-2} 



X uu ({2};f m _ 2 ,f m _ 2 ;{2iy,iy n - 2 }) 



{f(ab)}{f(ab)} 



0(0 + 1) 



+(3 - 2*/) 



(0 + 2) 



{/(aa)}{/(aa)} 
{/(aa)}{/(66)} 

{f(aa)}{f(ab)} 



o 



A1 « li b 



a) 6 "' 2 + 



(o + i)(o + 2) 



5^,2 



(r afe -l) 2 1 + (r a6 + l) 2 1 



4, 



(0 + 2) 



2(0 + a i) 



r„ b (r a6 + 1) n[ a) ^ab(r a6 - 1) u a b) ^ 

2(0 + 1) (0 + l)(0 + 2). 2(0 + 1)' 

-<V,2 ~ <Vi 



(0 + 1) 



2n- 
o 

(0 + 2) 







V,2 







(0 + 2) 



^ + ( 3-2,) ^ + 1 )^I 1 ) -^±I)^ 1 



(r ab - i)nj b) 







{/m-2} {/m- 2 } 



n/c/(/ m -2,/^ 2 ;{2,i^ 2 }) 



(O 2 - 1)1 1/2 



{/(a6)}{/(a&)} 
{/(a6)}{/(ac)} 
{/(a6)}{/(aa)} 




2 



(0 2 -4)_ 
O 



1 + Wni 6 ) 



+ 1- — 



(o 2 - 1) 



1/2 



2 L(0 2 -4) 

-o l ^-V 

[(02 -4) J 







T acJ 


n£ 6) 




1/2 J x 






~ 







O 



Finally, for f m -2 = {m — r — l,r — 1} irrep [this corresponds to (ii) in Fig. 3] we have 



r ab = m - 2r + 1 = 2F + 1 , 

r ai = m - r + i - 1, r bi = r + i - 2 ; i = 3, 4, . . . , Q, 

^ = (m-r + 1) (33) 

(m — r + Q — 1) 
r (o) = (r) 
b (r + Q-2) 
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These and A// m _ 2 /A/ m will give the formulas for the lower order moments of one and two 
point functions as described in Section 2. The dimension ratios needed are, 

M{m-r-2,r} (m — r) (m — T + 1 ) (iTL — 2f — 1) 



{» 



-■I 



m(m — l)(m — 2r + 1) 



{m— ) — l,i — 1} 
A/" [m— r,r} 



r(m — r + 1) 
m(m — 1) 



(34) 



M{ m - r ,r-2} r(r- l)(m-2r + 3) 



•N"{m-r,r} m(m — l)(m — 2r + 1) 

Using Eqs. (31)-(34) and the expressions in Table I, it is possible to derive analytical formulas 
for the P's, Q's and 7£'s that define (P 2 ), En and £22- The final formulas (obtained using 
MATHEMATICA) are, with (m, F) defining / m , 

P^ 2 >(m, F) = -- [3m(m - 2) + 4F(F + 1)] , 
P^ l2 >(m, F) = -- [m(m + 2) - 4F(F + 1)] , 
<2"=°({2}:m,F) = [P^(m,F)] 2 , 



p{ l2 >(m,F) 

(n + i) 



Q^=°({l 2 }:m,P) 
^m^F). i6(fl + 2) 

x [2(fl - 2)p{ 2 >(m, P) {3(2fi + m)(m - 2) + 4P(P + 1)} + 80(m - + 2m - 4)P(P + 1)] , 

(0 - l)P {l2} (m F) 
^^({l 2 } : m, P) = ^ '- K —^-t [{20. + m)(m + 2) - 4P(P + 1)] , 



Q^= 2 ({2}:m,P) 



(n) 



8(0 + 2) 



[(30 2 + 70 + Q)[F(F + l)f 



+ — m{m - 2)(20 + m)[20 + m + 2){0 -l)(0- 2) 

+ F{F+Y)_ + m ^ 5Q + _ ^ + _ ^ n _ 



Q^= 2 ({1 2 } :m,F) 



0(0-3)P^(m,F) 
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[(20 + m)(20 + m - 2) - 4P(P + 1)] , 



m, P) = P {2} (m, P) P^ l2 >(m, P) , 



n v =\m,F) 



' 2 -l (2-Q)P^(m,F) 
2 -A 8 



{4[P(P + 1) - 30] + 3m(20 + m - 2)} . 



(35) 



Note that Eq. (35) is related to the EGUE(2)-SC/(2) results given in [6] by the 0^-0 
transformation and they are also closely related to the results for spectral variances given 
in [8]. 
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(1) 



(ii) 



(iii) 



f m ={m-r,r} U(ty irreps 
— m-r - 



a a 



fm-2 ={m-r-2,r} 

















• 








• 


a 



fm-2 ={m-r-l,r-l} 



b b 

L-2 ={m-r,r-2} 



FIG. 3. Young tableaux denoting the two-rowed SU{Vt) irreps f m = {m — r, r} appropriate for 
BEGUE(2)-S'f/(2). Removal of two boxes generating m — 2 particle irreps f m -2 are also shown in 
the figure. For (ii) both the irreps f 2 = {2} and {l 2 } will apply while for (ii) and (iii) only {2} 
will apply. 

V. EMBEDDED GAUSSIAN UNITARY ENSEMBLE FOR SPIN ONE BOSONS: 
BEGUE(2)-SU(3) WITH r = 3 

Spin one boson systems, discussed in [9], possess U(3Q) D U(Q) <S> [SU(3) D SO(3)] 
symmetry. Instead of BEGOE(2) or BEGUE(2) generated by random two-body interactions 
preserving total spin S, it is also possible to consider interactions preserving the SU(3) 
symmetry. Then, for the GUE version, we have BEGUE(2)-5'?7 (3) that corresponds to r = 3 
in Section 2. As U(3) irreps will have, in young tableaux representation, maximum 3 rows, 
the U(Q) irrep also will have maximum three rows. Given m bosons in Q number of sp levels, 
the allowed U (ft) irreps are {/ x , f 2 , f 3 , / 4 , . . . , f n } with f x + f 2 + f 3 = m, f x > f 2 > f 3 > 
and fi = for i — 4, 5, . . . , ft. Because of the last condition we use simply {/i, f 2 , /s}. For 
f 2 — and / 3 = 0, we have totally symmetric irreps with {fi} = {m} and then all the results 
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(i) 


U(Q) in 

^- — r - — > 


eps 




















































t 


• 


a a 

fm ={r,r,r},f m _ 2 ={r,r,r-2}, f 2 ={2} 
(ii) 




















b 
a 
















• 
















t 


/ 


n ={r,r,r}J m -2 ={r,r-l,r-l}, f 2 ={l 2 } 



FIG. 4. Young tableaux denoting the three-column SU(Tl) irreps f m = {r, r, r}, m = 3r appropriate 
for BEGUE(2)-S'L r (2). Removal of two boxes generating m — 2 particle irreps f m -2 are also shown 
in the figure. For (i) only the irrep f 2 = {2} will apply while for (ii) only {l 2 } will apply. 

derived in Section 3 will apply directly. Similarly, for /2 ^ and f 3 = 0, all the results 
of Section 4 will apply. Thus the non-trivial irreps for BEGUE(2)-S77(3) are the m-boson 
irreps f m = {fi, f 2 , f 3 } with f 3 ^ 0. Given a f m in general there will be six f m - 2 and they 
are {/ a - 2, f 2 , f 3 }, {f u f 2 - 2, / 3 }, {f u f 2 , f 3 - 2}, {f, - 1, f 2 - 1, / 3 }, {f, - 1, f 2 , f 3 - 1}, 
{/i, f 2 — 1, f 3 — 1}. Therefore, as seen from Section 2, deriving analytical formulas for P's, 
Q's and TZ's that determine (H 2 ), En and E22 will be cumbersome. One situation that is 
amenable to analytical treatment is for the irreps {r, r, r}, m = 3r. For this class of irreps, 
the f m -2 are simple as shown in Fig. 4. For f m - 2 = {r,r,r — 2} we need lT a and II" and 
they are given by, 

it = 3r n" = 6r(r - 1} m) 

a fi + r-3' a (fi + r-3)(fi + r-4) ' 1 ' 

Similarly, for f m -2 — {r, r — 1, r — 1} we need r a b, 11^ and 11^ and they are, 

r .__i n w - 3r n (a) - 2(r + 1) f37) 

Tab " X ' iia _ 2(n + r-3)' 6 "(O + r-2)- 1 J 
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□ 
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<C 2 [SU(3)]> 



FIG. 5. (a) Variation of spectral widths as a function of m with fixed f m . Shown are the results 

(k) Ik) 

for f m as one-rowed f m irreps (red circles), two-rowed f m irreps (blue squares) and three-rowed 
fm^ irreps (magenta triangles), (b) Variation of spectral widths as a function of f m with fixed 
m. Shown are the results for $7 = 6 and m = 10, 15. Instead of showing f m , we have used 
(C 2 [SU(3)}) f ~"\ 



In addition, ratio of the Sq dimensions needed are, 



■A/"r,r,r-2 _ 2(r - 1) W r ,r-l,r-l __ T + 1 

M r ,r,r ~~ (3r - 1) ' Af r ,r,r _ (3r - 1) 
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(38) 



With these, carrying out simplification of the formulas given in Table I will give the following 
results, 



p{ 2 }(m, {r, r, r}) = — 3r(r — 1) , P^ l2 ^(m, {r, r, r}) = — —r(r + 1) , 

Q U=Q {{2} : m, {r, r, r}) = (3r) 2 (r - l) 2 , 

(3r) 2 (r+ l) 2 



Q^°({1 2 } : m, {r,r,r}) = 



Q u= \{2} : m,{r,r,r}) = 
Q^ =1 ({1 2 } :m, {r,r,r}) 



Q I '= 2 ({2} : m,{r,r,r}) 
Q^= 2 ({l 2 }:m, {r,r,r}) = 



6(tt + l)(n- 3)r(r- l) 2 (^ + r) 

(0+2) ' 
3(0 - l)(fl-3)r(r + l) 2 (^ + r) 

" 2(Q — 2) ' 

3fi(fi - 2)(fi - 3)r(r - + r)(fi + r + 1) 

4(fi + 2) 

3fi(fi - 3)(Q - 4)r(r + + r)(Q + r - 1) 



1Z U °(m, {r, r, r}) = 



(3r) 2 (r 2 - 1) 



8(Q-2) 



7^ (m, {r, r, r}) = — -t 



'Q 2 -l 

n 2 -4 



3(fi-3)r(r 2 - l)(fi + r) 



(39) 



Using these equations one can calculate the variances (H 2 ) and the covariances En and £22 
for irreps of the type {r, r, r}. For example, Eq. (15) can be simplified using Eq. (39) to 
give a compact formula for spectral variances, 



(H 2 ) m ' {r > r ' r} = A 2 2} 



+ A 2 l2} 



-r(r - l)(n + r - 3)(n + r - A) 



-r(r + + r - 2)(fi + r - 3) 



(40) 



It is also possible to derive analytical results for the irreps f m = {r+1, r, r} and {r+2, r, r} 
just as it was done for {4 r ,p} irreps for EGUE(2)-S77(4) ensemble in [7]. The results are as 
follows. For these irreps, the allowed f m -2 irreps and the corresponding r and IT ., functions 
are given in Table II and the dimension ratios in Table III. Using these, for f m = {r + 1, r, r} 
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FIG. 6. Self and cross correlations in energy centroids and spectral variances as a function 
of to and m' for £7 = 6 (with fixed f m and f m ')'- ( a ) [Sn(m, f m ; m', fm')] 1 / 2 with {/ m } = 
{to/3, to/3, to/3} for m mod 3 = 0, {f m } = {(to + 2)/3, (m — l)/3, (m — l)/3} for m mod 3 = 1 
and {fm} = {(rn + 4)/3, (m — 2)/3, (m — 2)/3} for m mod 3 = 2 and similarly / m / is de- 
fined; (b) [£ 22 (m, /mSm',/^)] 1 / 2 with {f m } = {m/3, m/3, to/3} for m mod 3 = 0, {/ m } = 
{(m + 2)/3, (m-l)/3, (m-l)/3} for m mod 3 = 1 and {f m } = {(m + 4)/3, (to - 2)/3, (m - 2)/3} 
for m mod 3 = 2 and similarly f m > is defined; (c) [Sn(m, f m ; m', fm')] 1 / 2 with {/ m } = {m/2, m/2} 
for m mod 2 = and {f m } = {(m + l)/2, (m — l)/2} for m mod 2 = 1 and similarly f m i is 
defined; (d) [E 22 (m, / m ; m', / m ')] 1/2 with {/ m } = {m/2, m/2} for m mod 2 = and {/ m } = 
{(m + l)/2, (m — l)/2} for m mod 2 = 1 and similarly f m i is defined; (e) [En(m, f m ; to', fm 1 )] 1 / 2 
with {f m } = {m} and {/ m /} = {to'}; (f) [£ 22 (m, / m ; to', / m ')] 1/2 with {/ m } = {to} and 

{/ TO '} = {"»'}• 
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FIG. 7. Self and cross correlations in energy centroids and spectral variances as a function of f m 
and f m i for Q = 6 (with fixed m = m'). Results are shown for (a) m = m' = 10 with f m = {10} 
(red circles), {5,5} (red stars) and {4,3,3} (red squares) with all one, two and three rowed f m i 
irreps; (b) m = m' = 15 with f m = {15} (blue circles), {8,7} (blue stars) and {5,5,5} (blue 
squares) with all one, two and three rowed f m i irreps. 
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irreps, P, Q and 1Z functions are, 

pW{m, {r + 1, r, r}) = -r(3r - 1) , P^fm, {r + 1, r, r}) = ~(5 + 3r) 
Q»=°{{2} :m,{r + 1, r, r}) = r 2 (3r - l) 2 , 

2 

Q^ =0 ({1 2 } :m,{r + 1, r, r}) = — (5 + 3r) 2 , 



rfl + ft") 

Q u=1 ({2} :m,{r+ 1, r, r}) = K - - - J [6r 3 (-3 + ft) + 3(-3 + ft)ft 
+2r 2 (-3 + ft)(-2 + 3ft) + r {-2 + 3(9 - 2ft)ft}] , 
Q^ =1 ({1 2 } : m, {r + 1, r, r}) = ^ + ^ [-r(5 + 3r) 2 
+ {-18 + r(-18 + r + 3r 2 )} ft + 3(1 + r)(2 + r)ft 2 ] , 

e-({2} : m , { r + 1, , r}) = ^ ±|^± + g> [8 + 3^-2 + «) <«> 
-(-8 + ft)ft + r(-2 + ft)(l + 3ft)] , 

Q- 2 ({1 2 } : m, {r + 1, r, r}) = r( ~ 3 ^2+^ ~ + ~ l~ 16 + ^ A + ^ 
+r(-4 + ft)(7 + 3ft) + ft(-14 + 5ft)] , 
n u=0 {m, {r + l,r,r}) = — (-1 + 3r)(5 + 3r) , 

^ =1 (m, {r + 1, r, r}) = i r ^ m ~ 151Q ^ + 153 (~ 3 + ^ Q 

+r 2 {600 + (283 - 151ft)ft} - 5r {60 + 13ft(-7 + 2ft)}] . 

Similarly, for / m = {r + 2, r, r} irreps we have, 

p{ 2 >(m, {r + 2, r, r}) = -(3r 2 + r + 1) , p{ l2 }(m, {r + 2, r, r}) = --(7 + 3r) , 

Q^ =0 ({2} : m, {r + 2, r, r}) = (3r 2 + r + l) 2 , 

Q-°({l 2 }:m,{r + 2,r,r}) = ^(7 + 3r) 2 , 

Q" =1 ({2} : m, {r + 2, r, r}) = [-4(1 + r + 3r 2 ) 2 

+ {2 + r {3 + r {51 + 4r(-7 + 3r)}}} ft + (2 + llr + 12r 3 )ft 2 ] , 

Q^({1 2 } : m, {r + 2, r, r}) = ± | [-r(7 + 3r) 2 

+ (5 + 3r)(-6 + r 2 )ft + {10 + 3r(4 + r)} ft 2 ] , 

Q- 2 ({2} : m, {r + 2, r, r}) = " [3r 4 (-3 + ft)(-2 + ft) + (-1 + ft)ft(2 + ft)(3 + ft) 

+2r 3 (-3 + ft)(-2 + ft)(4 + 3ft) + r 2 (-3 + ft)ft {7 + 3ft(l + ft)} 
+r(-3 + ft) {22 + ft {42 + ft(19 + ft)}}] , 

(42) 
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TABLE II. Formulas for the functions H^__ ^'s denned in Eq. (23) as required for {f m } irreps 

{r + l,r, r} and {r + 2,r, r}. Given also are the values of axial distances (r 's). Also, 11^ = 

r/(Q + r — 3) for both {f m } examples shown in the Table. 
{fm} {fm-2} Required functions 

{r + l,r,r} {r + l,r,r-2} Jl> a = , IL>> - 5 '' ( '' " ' 



3(fi + r-3)' a (ft + r-3)(0 + r-4) 
/(aa) 

\r \ \ r It 1} r h - 1 II (fe) - ^ I4 (a) - 3 ( r + T ^~ 

/(a&) 

r 11 o rr( c ) 2r rr(a) ^ + 3 

f(ac) 

{r + 2 ,r,r} {,■■■>.,:,■ 2} 1% = — - , Iff - 9 '' { '' " " 



2(0 + r- 3) a 2(n + r -3)(fi + r -4) 
/(aa) 

r r 1 2 r 1 r llr,- 1 II (fe) - 5r n (a) - 4( r + 1} 

(r + 2,r l,r 1} r ab - 1, ll„ - 4(n + r _ 3) > "ft " 3(0 + r - 2) 

/(at) 

/ r + i ~ r _u - _ 4 n (c) - 2r n (a) - 2(r + 4) 

Hoc) 

~ " ~. r + 4 ~ (r + 3)(r + 4) 

n e = TT7?n r-TT' n c- 



2(fi + r + l)' c 6(0 + r)(0 + r + l) 



/(cc) 



Q- 2 ({1 2 } : m, {r + 2, r, r}) = r( 3 + ^ + ^ - + ^ [~ 44 + 3r 2 (-4 + fl) 

+r(-4 + + 3fi) + Q(-12 + 7Q)\ , 

<Jl»=°(m,{r + 2,r,r}) = Ul + 3r)(l + r + 3r 2 ) , 

^ =1 (m, {r + 2, r, r}) = [8(7 + 3r)(12 + r + 37r 2 ) 

+ {-1008 + r {1106 + (241 - 289r)r}} Q + {208 - r(489 + 289r)} Q 2 } . 

In addition to analytical results, as stated in Section 2, one can use the tables in [7] (also 
to a large extent Table I) and obtain numerical results for the variation of spectral variances 
with the eigenvalues of the quadratic Casimir invariant of U(Q) or equivalently C 2 [SU(3)}, 
for various (Q, m) values and also for both self and cross correlations in energy centroids 
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TABLE III. Dimension ratios with respect to the S m group for the examples in Table II. 

{fm} {fm-2} }r m ^ 

~r 7~r i 5(r - 1) 

{r + l,r,r} {r + l,r,r-2} - ± 

< r + 1 ' r - 1 ' r - 1 >3k + ^ 
/ u r + 3 ) 

{r ' r ' r ~ 1} gE+n 

{r + ^r.r} {r + 2,r,r-2} 2( ^ ( ^ J + 1) 

( r + 2 r -l r-U 5r (r + l) 

i + ' ' 1 3(3r + 2)(3r + l) 

4r(r + 4) 

{r + 1 ' r ' r " 1} 3(3r + 2)(3r + l) 
r ; ( r + 3 )(r + 4) 

ir ' r ' r| 6(3r + 2)(3r + l) 

and spectral variances. Note that For a Q = 6 system with A? 2 i = A{i2} = 1 calculations are 
carried out for various choices of m and f m and the results are shown in 5, 6 and 7. Let us 
mention that C 2 [SU(3)] for a irrep {f\, f 2 , fs} is given by the formula, 

(C 2 [SU(3)]) {fM} = A 2 + /i 2 + A/i + 3(A + fi) ; A = A - / 2 , = / 2 - / 3 . (43) 

It is seen from Fig. 5a that the spectral widths will be largest for one rowed irreps and 
smallest for three row irreps for a fixed m. Also, widths as expected increase with m. 
Similarly, Fig. 5b shows that for a fixed m, widths increase as the eigenvalue of C2[SU(3)] 
increases and this is consistent with the observation in Fig. 5a as the eigenvalue of C 2 [SU (3)] 
is largest for totally symmetric irrep. Results in Fig. 6 show that: (i) the centroid and 
variance fluctuations increase with w! for fixed m and vice-versa; (ii) they are larger for 
three rowed irreps compared to those for one rowed irreps; (iii) centroid fluctuations are 
much larger than variance fluctuations as seen before also for EGUE(2)-s and EGUE(2)- 
SU(4) ensembles. Similar trends are also seen for m — m! but varying f m i with fixed f m and 
these results are shown in Fig. 7. A different trend is seen for the covariances in spectral 
variances for the totally symmetric irrep f m = {m} and f m > varying. More importantly, the 
centroid and variance fluctuations are smallest for the ground state i.e., the most symmetric 
irrep for bosons. It is seen from Figs. 6 and 7 that the covariances in energy centroids are 
~ 15 — 25% and the covariances in spectral variances are ~ 8 — 15%. 
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VI. CONCLUSIONS 



In this paper, given first is a general formulation for deriving lower order moments of 
the one- and two-point correlation functions in eigenvalues that is valid for any embedded 
random matrix for fermions as well as for bosons with U (Q) (g> SU (r) embedding and with 
two-body interactions preserving SU (r) symmetry. Results of the present paper unify all the 
results known before for EGUE(2)'s and BEGUE(2)'s. Presented are new results for boson 
systems with SU(r) symmetry, r = 2,3. These results should be useful in future studies 
of two species boson systems and spin one boson systems. In future, it will be useful to 
derive analytical forms for SU (Q) Racah coefficients [6, 7] or develop tractable methods for 
their numerical evaluation to establish Gaussian form of the eigenvalue densities generated 
by embedded ensembles with SU(r) symmetry both for boson and fermion systems. 
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